Time dependence of laser cooling in optical lattices 
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We study the dynamics of the coohng of a gas of caesium atoms in an optical lattice, both 
experimentally and with ID full-quantum Monte Carlo simulations. We find that, contrary to the 
standard interpretation of the Sisyphus model, the cooling process does not work by a continuous 
l/^ , decrease of the average kinetic energy of the atoms in the lattice. Instead, we show that the 

■ momentum of the atoms follows a bimodal distribution, the atoms being gradually transferred from 

' a hot to a cold mode. We suggest that the cooling mechanism should be depicted in terms of a rate 

, model, describing the transfer between the two modes along with the processes occurring within 

^ . each mode. 



X 



PACS numbers: 32.80.Pj, 05.10.Ln 



I. INTRODUCTION 



Laser cooling is now a well-established technique enabling to reach temperatures of the order of the microkelvin in a 
Q ' gas of atoms Q. It has made possible Bose- Einstein condensation 0, while laser-cooled atoms are used, for instance, 
in atomic clocks Q and as frequency standards Q , and have been proposed for parity violation measurements and 
quantum information processing . 

An appropriate superposition of laser beams can result in a spatially-periodic modulation of the polarisation of the 
light, creating an optical lattice 7, 8, 9]. A multilevel atom with a degenerate ground state moving through such a 
lattice will experience large non-adiabatic effects as the light shift of the different substates and the optical pumping 
rates then depend on the position of the atom. The standard model used to explain laser cooling in optical lattices 
rS i' is Sisyphus cooling J.Q., 11], where an atom will preferentially jump from the crest of a potential hill to the valley 
of an other potential, losing kinetic energy each time. Atom dynamics can be expressed in a semi-classical way by 
■ this model as the interplay between a velocity-dependent cooling force due to the Sisyphus effect and a diffusion term 
^ ' describing the recoil from photon absorption-emission cycles along with fluctuations in the gradient force when going 
- from one potential curve to the other [lOl [ll| . 

Although it was predicted many years ago , it is only recently that experimental observations have shown non- 
Gaussian velocity distributions in optical lattices iTp , as well as in atomic fountains . While the standard Sisyphus 
model considers a cooling force linear in velocity 11], leading to a Gaussian momentum distribution, refinements using 
non-linear forces give non-Gaussian distributions, following for instance Tsallis or Levy statistics [15.. 16j. 

Results obtained by taking into account the localisation of the atoms around the potential minirna also show deviations 
from a Gaussian |17lll8l |. Moreover, arguments based on a "band" picture of optical lattices 1!||, where some atoms 
^ are trapped in the potential wells of the lattice while others move around in the "conduction band" , call for a bimodal 
• I— I ■ distribution [13) Indeed, it was found that a double Gaussian function provides the best fit to both experimental 
and simulated data for the steady state momentum distribution of atoms [l^ . 
(— I ' We report in this letter a study of the time evolution of the cooling of a gas of caesium atoms inside an optical 
^ y lattice. It provides new elements to understand the origin of the non-Gaussian momentum distributions. Apart from 
' a few studies concerning localisation 0,1121, diffusion and damping rates [3 of atoms in optical lattices, 
this is one of the first observation of the dynamics of the cooling process. 



II. METHODOLOGY 



Details of the experimental setup can be found in refs. |2(i l27l |. Briefly, we first accumulate ^'^^Cs atoms in a 
magneto-optical trap. We adjust the irradiance and the detuning, then we turn off the magnetic field and leave the 
atoms in an optical molasses with even further reduced irradiance and increased detuning. The atoms are thus cooled 
to about 40 /xK before being transfered to the optical lattice; such a "hot" initial cloud of atoms is used to reduce the 
effects from the initial kick given to the atoms when the lattice is turned on. The lattice itself has a 3D structure 
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made up from a four-beam configuration [7], |9j : two laser beams are linearly polarised along the x axis and propagate 
in the yz plane symmetrically with respect to the z axis, whereas the other two beams are polarised along the y 
axis and propagate in the xz plane symmetrically with respect to z. This yields a tetragonal pattern of points with 
pure circular polarisation, alternately cr+ and and a~ . The lasers are slightly detuned below the resonance of the 
Fg = 4 — > Fc = 5 transition of the D2 line (6s^Si/2 ^ 6p^P3/2) at A = 852 nm. After a given time the optical lattice 
is turned off (in 4.2 /is) and the atoms are left to fall and then detected by a probe beam located ~ 5 cm below the 
trapping region. This turn-off time is fast enough so that no adiabatic cooling of the atomic cloud can occur. The 
momentum distribution of the atoms is recovered from the time-of-flight signal. 

The simulation of the dynamics of the atoms in the optical lattice is carried out using a full-quantum Monte Carlo 
wave- function method for the = A ^ = h transition of caesium • Due to the numerical complexity of the 
3D case, we restrict ourselves to a ID lin _L lin configuration, reproducing the same alternation of and a~ potential 
wells (see, e.g., ref. Q). The time dependence of the cooling process is obtained by considering a series of different 
histories where the initial momentum of the atom is chosen randomly from a normal distribution corresponding to a 
temperature of 50 ^K. Observables are recovered by averaging over these histories, e.g., the momentum distribution 
is obtained from the time-dependent wave functions iph by 

1 ^ 

D{p:t) ^ -Y.\'i^H{p;t)\^ , (1) 

h=l 

where the index h labels the different histories. 



III. RESULTS 



Sample momentum distributions obtained by the numerical simulation are given in fig.H for a detuning of A = — lOF 
and a potential well depth of [/o ~ 127 E^, with T/2tt = 5.2227 MHz the natural linewidth of the 6p^P3/2 level of 
caesium and E^ = huj^ the recoil energy, where cji./27r — 2.0663 kHz is the atomic recoil frequency '2^. The momentum 
is expressed in units of the recoil momentum — h/X. Starting from the initial distribution at 50 /iIC, we clearly see 
that a narrow peak gradually grows in the centre of the momentum distribution. At the same time, the population 
of the wings from the initial distribution decreases and spreads out to higher momentum. It appears that the atoms 
are not gradually cooled but populate progressively the central feature: they are transfered from a "hot" to a "cold" 
mode. To analyse these results, we have fitted the momentum distribution to a double Gaussian, 

Dip) - v4eoid exp [V/ (2a2„id)] + ^hot cxp [-p^ (2aLt)] , (2) 

by finding the parameters (A and cr) that best reproduce the second, {p'^}, and fourth, (p^), moments of the dis- 
tribution while conserving the norm, with the constraint CTcoU < o'^ot- This method appears better suited than the 
traditional fitting procedure to reproduce the wings of the distribution when the cold mode dominates. We find, 
fig.[2Ia), that there is no time evolution of the width acow, hence we can assign a temperature Tcom — '^coid/ i''^^^) 
to the cold mode, where m is the mass of the atom and fee the Boltzmann constant. This is another striking feature 
of the cooling process, which works not only by transferring atoms between the two modes, but puts them in an 
unchanging cold mode. On the other hand, atoms in the hot mode are continuously heated. We have checked that 
no steady state is obtained in the hot mode for longer times. For as long as we can see, the momentum distribution 
in the hot mode shows a Gaussian profile in the region where it can clearly be distinguished from the cold mode. 
In contrast, as can be seen in fig. |3b), the relative population of each mode (calculated as \/2TTAa for a normalised 
distribution) reaches a steady state after some time. By fitting the population of the cold mode to an exponential 
function, we recover a transfer rate of « 7.9 x 10^ s~^. This rate is about 15 times slower than the localisation 
rate previously measured for caesium in a ID lattice [21} and 2 orders of magnitude slower than the scattering rate. 
The persistence of hot atoms seems to be due to the fact that some higher momentum atoms are heated and never 
cooled. Indeed, starting from a lower temperature such as 5 /iK leads to a greater transfer from the hot to the cold 
mode, with less heating in the hot mode. 

Turning now to the experimental results, we show in fig.|3|the momentum distribution obtained for a laser detuning 
of A = — 12.6r and an irradiance corresponding to a potential depth Uq — 217 E^. (Since the simulation is for a ID 
lattice, a quantitative comparison to the experiment is not meaningful, so no attempt is made to exactly match the 
parameters for both cases.) We see again the appearance of a central feature of almost constant width. The main 
difference with the simulation is that after a certain time the signal decreases, as seen in fig. Elb). This is due to an 
overall loss of atoms due to their diffusion out of the lattice and collisions with background gas, while the numerical 
simulation considers an effectively infinite lattice. This phenomenon is probably also responsible for another difference 
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between the two results, namely the absence in the experimental signal of a significant heating of the hot mode, as 
the high momentum atoms seen in the simulation would quickly escape from the lattice region. 

The result obtained by fitting the experimental data to a double Gaussian is given in fig. 01 Where have reverted 
here to a method as a fit using the moments of the distribution gave aberrant results due to the bump present in 
the experimental data at negative momentum. It is clear from the time evolution of Chot (although it is noisy) that 
there is no heating observed in the hot mode. This is nevertheless not necessarily contradictory with the result of 
the simulation, as the decrease ahot could be explained by the preferential loss of the high momentum atoms. The 
temperature in the cold mode is almost constant, the slight decrease with time being within the experimental and 
fitting errors. The relative populations of the two modes shown in fig. ^b) seem to indicate that a steady state is 
not reached, with the cold mode continuously gaining in population, which would indicate that the loss of atoms is 
through the hot mode. The rate of transfer to the cold mode is found to be of the same order as for the simulation, 
with « 5.8 X 10^ s~^, which is 6 times slower than the localisation rate in a 3D lattice '2l'| and ~ 400 times slower 
than the scattering rate. 

We ran simulations over a range of the parameters A and Uq and the temperatures thus obtained for the cold 
mode are given in fig^El The variation of TcoM as a function of the well depth looks similar to that obtained in 
previous studies |l2l ll£L .23] , although the temperature was then calculated from the root-mean-square momentum. 
It also appears that the temperature does not follow the linear dependence with potential depth predicted by an 
analytical model of Sisyphus cooling 's^l for the entire range considered here. We also note that no decrochage is 
seen, i.e., there is no threshold potential below which there would be a rapid increase of the temperature. The curve 
of fig.|3a) is in fact similar to the results previously obtained by considering not the r.m.s. momentum, but the width 
of the distribution at 1/ -^e ■ The dominance of the cold mode over the hot mode is such that the width at 

l/^/e is essentially equal to (Tcoid, except for very short times. It thus seems that the previously observed decrochage 
appears when the population of the hot mode is great enough to influence significantly the r.m.s. momentum. It is a 
feature of the relative population of the modes, not of the cooling temperature reached. Extrapolating to a vanishing 
potential depth results in a non-zero value of Tcoid, but we expect the behaviour to become non linear as the number 
of bound states in the lattice goes to zero. The temperature of the cold mode also varies with the detuning, as shown 
in fig. EJb). This is consistent with what has been observed for rubidium |23l l25l|. although it appears that Tcom 
decreases exponentially with |A|, instead of following the functional form proposed in refs. pa. l25|. 

IV. DISCUSSION 

Our results call for a discussion of the cooling of atoms in optical lattices in terms of rates for a bimodal distribution, 
based on the hypothesis that the cold and hot modes correspond to trapped and untrapped atoms, respectively. The 
dominant process during the initial cooling stage is the transfer from the hot to the cold mode. Within the hot mode, 
simulations show that there is an important heating effect, while for the cold mode basically no heating or cooling is 
seen, which means that there is a fast equilibrium reached inside the cold mode. This latter process depends on the 
detuning of the lasers, as seen by the dependence of T^oid not only on the potential depth but also on the detuning. 
There is also the transfer from the cold to the hot mode, which, according to simulations, is slow with respect to 
the transfer from hot to cold. Otherwise, there would be a gradual depletion of the cold mode as atoms in the hot 
mode are heated to the extent where they cannot go back to the cold mode. Nevertheless, under the hypothesis 
that cold atoms are trapped, the loss of atoms due to diffusion could only be explained by a coupling from cold to 
hot. We also note that associating the cold mode with trapped atoms is not incompatible with the observation that 
localisation is faster than cooling, as atoms can spend more time near the potential minima {i.e., localise) while still 
being untrapped. In addition, the rate actually measured in ref. [2ll | may also reflect the progressive filling of lattice 
sites. 



V. CONCLUSION 

To summarise, we have shown that atoms in an optical lattice follow a bimodal distribution of hot and cold atoms. 
The cooling inside the lattice takes place via a transfer of the atoms from the hot to the cold mode. We suggest 
that the cooling of atoms in an optical lattice be depicted in terms of a rate model, describing the transfer between 
the two modes along with the processes occurring within each mode. This calls for further investigations in order to 
understand by exactly which mechanism this is taking place. 



4 



Acknowledgments 

We thank R. Kaiser, E. Lutz, K. M0lmer, and L. Sanchez-Palencia for stimulating discussions. This work was 

supported by Carl Tryggcrs stiftclsc, Magnus Bcrgvalls stiftclsc, Kcmpcstiftclscrna, K & A Wallenbergs stiftelse, 
and the Swedish Research Council. S.J.H.P. thanks Carl Tryggers stiftelse for financial support. This research was 
conducted using the resources of the High Performance Computing Center North (HPC2N). 



[1] S. Chu, Rev. Mod. Phys. 70, 685 (1998); W. D. Phillips, ibid. 70, 707 (1998); C. Cohen-Tannoudji, ibid. 70, 721 (1998). 

[2] E. A. Cornell and C. E. Wieman, Rev. Mod. Phys. 74, 875 (2002); W. Kettcrlc, ibid. 74, 875 (2002). 

[3] S. Bize, P. Laurent, M. Abgrall, H. Marion, I. Maksimovic, L. Cacciapuoti, J. Griinert, C. Vian, F. Pereira dos Santos, P. 

Rosenbuscli, P. Lemonde, G. Santarelli, P. Wolf, A. Clairon, A. Luiten, M. Tobar, and C. Salomon, J. Phys. B: At., Mol. 

Opt. Phys. 38, S449 (2005). 

[4] L. HoUberg, C. W. Gates, G. Wilpers, C. W. Hoyt, Z. W. Barber, S. A. Diddams, W. H. Oskay, and J. C. Bergquist, J. 

Phys. B: At., Mol. Opt. Phys. 38, S649 (2005). 
[5] S. Sanguinctti, J. Gucna, M. Lintz, P. Jacquier, A. Wasan, and M.-A. Bouchiat, Eur. Phys. J. D 25, 3 (2003). 
[6] C. Monroe, Nature 416, 238 (2002). 

[7] P. Jessen and I. Deutsch, Adv. At. Mol. Opt. Phys. 37, 95 (1996). 
[8] L. Guidoni and P. Verkerk, J. Opt. B.: Quantum Semiclass. Opt. 1, R23 (1999). 
[9] G. Grynberg and G. Robilliard, Phys. Rep. 355, 335 (2001). 
[10] Y. Gastin, J. Dalibard, and C. Gohen-Tannoudji, in Light Induced Kinetic Effects on Atoms, Ions, and Molecules, edited 

by L. Moi, S. Gozzini, G. Gabbanini, E. Arimondo, and F. Strumia (ETS Editrice, Pisa, 1991), pp. 5-24. 
[11] C. Gohen-Tannoudji, in Les Houches, Session LIII, Systemes Fondamentaux en Optique Quantique/ Fundamental Systems 

in Quantum Optics, edited by J. Dalibard, J. M. Raimond, and J. Zinn- Justin (Elsevier, Amsterdam, 1992), pp. 1-164. 
[12] J. Jersblad, H. Ellmann, K. St0chkel, A. Kastbcrg, L. Sanchez-Palencia, and R. Kaiser, Phys. Rev. A 69, 013410 (2004). 
[13] Y. Sortais, S. Bize, C. Nicolas, A. Clairon, C. Salomon, and C. WiUiams, Phys. Rev. Lett. 85, 3117 (2000). 
[14] T. W. Hodapp, C. Gerz, G. Furtlehner, G. I. Westbrook, W. D. Phillips, and J. Dalibard, App. Phys. B: Lasers Opt. 60, 
135 (1995). 

[15] E. Lutz, Phys. Rev. A 67, 051402(R) (2003). 

[16] F. Bardou, J. -P. Bouchaud, A. Aspect, and C. Cohen-Tannoudji, Levy Statistics and Laser Cooling (Cambridge University 

Press, Cambridge, 2002). 
[17] J. Javanainen, Phys. Rev. A 46, 5819 (1992). 
[18] K. M0lmer and C. L Westbrook, Laser Phys. 4, 872 (1993). 
[19] Y. Gastin and J. Dalibard, Europhys. Lett. 14, 761 (1991). 

[20] L. Sanchez-Palencia, P. Horak, and G. Grynberg, Eur. Phys. J. D 18, 353 (2002). 

[21] G. Raithcl, G. Birkl, A. Kastberg, W. D. PhiUips, and S. L. Rolston, Phys. Rev. Lett. 78, 630 (1997). 
[22] G. Raithel, G. Birkl, W. D. Phillips, and S. L. Rolston, Phys. Rev. Lett. 78, 2928 (1997). 

[23] F.-R. Carminati, M. Schiavoni, L. Sanchez-Palencia, F. Renzoni, and G. Grynberg, Eur. Phys. J. D 17, 249 (2001). 
[24] F.-R. Carminati, L. Sanchez-Palencia, M. Schiavoni, F. Renzoni, and G. Grynberg, Phys. Rev. Lett. 90, 043901 (2003). 
[25] L. Sanchez-Palencia, M. Schiavoni, F.-R. Carminati, F. Renzoni, and G. Grynberg, J. Opt. Soc. Am. B 20, 925 (2003). 

[26] ,J. .Jersblad, H. Ellmann, and A. Kastbcrg, Phys. Rev. A 62, 051401(R) (2000). 
[27] H. Ellmann, J. Jersblad, and A. Kastbcrg, Eur. Phys. J. D 13, 379 (2001). 
[28] Y. Gastin and K. M0lmer, Phys. Rev. Lett. 74, 3772 (1995). 

[29] D. A. Steck, Cesium D line data. Technical Report No. LA-UR-03-7943, Los Alamos National Laboratory (2003), 

http : //steck . us/ alkalidata/. 
[30] J. Dalibard and C. Gohen-Tannoudji, J. Opt. Soc. Am. B 6, 2023 (1989). 



5 



3 



0.10 



0.08 



0.06 



0.04 



0.02 



0.00 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.0 ms 

0.1 ms 

; 0.5 ms 

: ; - - - - 1.0 ms 

5.0ms 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
80 -60 -40 -20 20 40 60 80 

P/Pr 



FIG. 1: Simulated momentum distribution [see eq. Q] at different lattice times, for A = — lOF and Uo = 127i?r. 




FIG. 2: Time dependence of the Gaussian fit parameters [see eq. (|5J] for the cold and hot modes, for the simulated data shown 
in fig. Q (a) widths CTcom and (Thot; (b) relative populations. 
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FIG. 3; Experimental momentum distribution at different lattice times, for A = — 12.6r and Uo — 217-Er- 
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FIG. 4: Time dependence of the Gaussian fit parameters [see eq. for the cold (squares) and hot (circles) modes, for the 
experimental data shown in fig. El (a) Widths CTcoM and ahot. (b) Relative populations, with fits to exponential functions (solid 
lines). Note that there is a scale change in the time axis. 
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FIG. 5: Temperature of the cold mode: (a) as a function of the potential depth, for detunings of —5T (squares), — lOF (triangles), 
and — 20r (circles); (b) as a function of detuning, for well depths of 127i?r (squares) and 196Ei (triangles). 



